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A semi- analytic approach to the computation of the non-linear power- spectrum of dark matter 
density fluctuations is proposed. The method is based on the Renormalization Group technique 
and can be applied to any underlying cosmological model. Our prediction on the baryonic acoustic 
oscillations in a ACDM model accurately fits the results of N-body simulations down to zero redshift, 
where perturbation theory fails. 



PACS numbers: 98.80.-k 

Cosmological linear perturbation theory is a funda- 
mental tool in the physics of the Cosmic Microwave Back- 
ground (CMB), allowing a careful extraction of cosmo- 
logical parameters from the data. In the case of Large- 
Scale Structure, the accuracy of theoretical predictions 
is not yet at the same level than for CMB, due to the 
higher degree of non-linearity of the underlying density 
fluctuations. For instance, the study of Baryonic Acous- 
tic Oscillations (BAO) - a fundamental probe for Dark 
Energy measurements - requires a precision of a few per- 
cent in the theoretical predictions for the matter power- 
spectrum (PS) in the wavenumber range k :^ 0.05 — 0.25 
h/Mpc [1]. Including higher orders in PT [2] is known 
to give a poor performance in this range, leaving N-body 
simulations as the only viable approach to the problem. 

Recently, however, perturbation theory (PT) has ex- 
perienced a renewed interest, mainly motivated by two 
reasons. First, next generation galaxy surveys are going 
to measure the PS at large redshift, where the fluctua- 
tions are still in the linear regime and 1-loop PT is ex- 
pected to work [1]. Second, Crocce and Scoccimarro [3] 
have shown that the perturbative expansion can be reor- 
ganized in a very convenient way, which allows the use of 
standard tools of field theory, like the Feynman diagrams. 
They managed to compute the two-point correlator be- 
tween density or velocity field fluctuations at different 
times (the 'propagator') by resumming an infinite class 
of diagrams at all orders in PT. Other approaches can be 
found in Ref. [4]. 

In this Letter^ we will proceed along the same path, by 
implementing Wilsonian Renormalization Group (RG) 
techniques to compute the PS. RG methods, widely used 
in statistical mechanics and quantum field theory [5] , are 
particularly suited to physical situations in which there 
is a separation between the scale where one is supposed 
to control the 'fundamental' theory and the scale were 
measurements are actually made. Starting from the fun- 
damental scale, the RG flow describes the gradual inclu- 
sion of fluctuations at scales closer and closer to the one 



relevant to measurements. The new fluctuations which 
are included at an intermediate step, feel an effective the- 
ory, which has been 'dressed' by the fluctuations already 
included. In the present case, the RG flow will start from 
small wavenumber s /c, where linear theory works, to reach 
higher and higher k. 

We will consider a self-gravitating system of Dark Mat- 
ter (DM) particles which, in the "single-stream" ap- 
proximation, is governed by the continuity, 9r ^ + V • 
[(1 + 5)m] = 0, and Euler equations dr v+7Yv+(v- V)v = 
— V^. Here r is the conformal time, 6 the mass-density 
fluctuation, v the peculiar velocity, (j) the peculiar gravi- 
tational potential which, on subhorizon scales, obeys the 
Poisson equation = 1 ^ , with H = dlna/dr, 

having assumed an Einstein-de Sitter background cos- 
mology. 

Following Ref. [3] we introduce the doublet (/?a(k, 77) 
(a = 1, 2) - defined in Fourier space - given by {cpi^ (^2) = 
{5^ — • v/7Y), where the time variable has been re- 
placed by the logarithm of the scale factor, 77 = ln(a/ain), 
ain being the scale factor at a conveniently remote epoch, 
when all the relevant scales are well in the linear regime. 
Notice that, compared with the definition in Ref. [3], we 
have an overall factor e~^, such that the linear growing 
mode corresponds to (pa = const. We can then rewrite 
the continuity and Euler equations in compact form, as 

{Sabdrj ^^ab)^b{Kv) = e'^7a6c(k, -p, -q) 

X ^6(p,^) V^c(q,^) , (1) 



where ft 



and repeated in- 



-3/2 3/2 J 

dices/momenta are summed/integrated over. 

We have defined a vertex function, 7a5c(k, p,q) 
{a^b^c^= 1,2) whose only non-vanishing elements are 
7121 (k, p, q) = 1/2 (^D(k+p+q)a(p,q), 7222 (k, p, q) = 
(5D(k + p + q)/?(p,q), and 7121 (k, p, q) = 7112 (k, q, p), 
where a(p,q) = [(p + q)-p]/p^ and /5(p, q) = (p + q)^(p- 
q)/(2p^g^). 
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Besides the vertex, the other relevant dynamical quan- 
tity is the linear retarded propagator, which gives the 
linear evolution in r] of cpa, ^aO^, rj) = gabiv^ vl^bO^^ 
with ri > ri' (the subscript "0" indicates solutions of the 
linear equations, e'^jabc = 0), which reads Qabiv^v') — 
[B + A exp(-5/2(7/ - vl)]abO{r] - i), with the step- 
function, B = I 3 2 ) and A = I ^ _^ ^ 

The growing {^pa oc const.) and the decaying {cpa oc 
exp(— 5/2/^)) modes can be selected by considering initial 
fields (fa proportional to Ua = (1, 1) and Va = (1, —3/2), 
respectively. 

To extend the validity of this approach to ACDM, we 
reinterpret the variable 77 as the logarithm of the linear 
growth factor of the growing mode, i.e. r] = \n{D~^ /D^). 
This approximation has been shown to accurately fit N- 
body simulations for different cosmologies [1]. 

Our aim is to apply methods familiar in quantum field 
theory to construct generating functionals for quantities 
like the PS, bispectrum, etc. ... The starting point is 
to write down an action giving the equation of motion 
(1) at its extrema. One can realize that a new, auxiliary, 
doublet field Xa has to be introduced to this aim, and 
that the action is given by 

S = J df] [Xa(-k, T]) {Sabdn + ^ab) ^b{K V) (2) 
-e'^Tabcl-k, -p, -q)xa(k, ri)Lpi,{p, r])ipc{q, r])] . 

Varying the action S w.r.t. Xa gives precisely Eq. (1), 
while varying it w.r.t. (pa yields an equation solved by 
Xa = 0. 

The probability of a classical field configuration is a 
delta function centered on the solution of the equation of 
motion. By averaging over Gaussian initial conditions, 
one can show [6] that the generating functional of corre- 
lation functions involving the cp and x fields is the path- 
integral 

Z[Ja, K5; P'] (3) 

VifaVxbexp^^j dr]dr]' [-\xaPabKv)Kv')Xb 

+ ^XaQab^h - ie"^ labcXa^b^c + iJa^a + ^^bXb] | , 

where Ja and are external sources and the momen- 
tum dependence is everywhere implicit. Here P^^^ik) = 
P^{k)uaUi) and P^ is the linear PS at the initial "time" 
?7 = 0, having assumed pure growing- mode initial condi- 
tions. As usual, we define also the generator of connected 
Green's functions W[Ja^Ki)] = —i\nZ[Ja,Ki)] and, for 
1-particle irreducible (IPI) ones, the effective action 
r[(^a, Xb] =W[Ja,Kij]- J dr]d^k{Ja(pa^Ki,Xb)/^^ terms 
of the "classical" fields, (paiJc^^d] = ^W[Jc, K(i]/6Ja and 
Xb[Jc^K(i] = SW[Jc^ Kd]/SKi,. In linear theory the path- 
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power spectrum: PatiVa, Vb; k) 



interaction vertex: —ie^jabcO^a: kb, kc 



FIG. 1: Feynman rules 



integral can be performed analytically, 
Zo[Ja,i^5;P'] = expj- J dr]dr]'[ya{v)Pab{v.v)Jb{v) 
+ iJa{v)9ab{v.v')Kb{v')]]. (4) 



with P^^{r]a^Vb;k) = gac{Va^O)gbd{Vb^O)P^^{k) the lin- 
early evolved PS. 

Standard methods of PT can be applied, using the 
Feynman rules sketched in Fig. 1. There are three build- 
ing blocks: the linear propagator gab^ the linear PS 
and the vertex e^^abc- The Feynman diagrams con- 
structed using these rules are in one-to-one correspon- 
dence with those of Ref. [3] , and reproduce all the known 
results of PT [2]. 

The full propagator, PS and vertex can be obtained 
from W and F as {6'^W/5JaSKi,)j^K=Q = -(^(k + 
kOGa6, {S^W/6Ja6J,)j^K=o = iS{k + k')Pab and 

{S^r/SXa^^b^^c)y:^, X=0 = -TabciK + k^ + T]). 

Using the definitions of effective action and classi- 
cal fields above, one can show that the full PS has 
the structure Pab = Pab + Pab ^ where Pl^{k]r],r]') = 
Gac{k',v^mM{k',i.^)PUk) and P,Y(fe;7?,V) = 
Q di' di"Gac{k', 7/, i')Gbd{k', V, i")^cd{k-, i"). 
where ^^ab is defined through {^'^T / 5xaXb)^,x=^ = 
iP'ab{mri)Ki)^^^ab- 

The starting point of our formulation of the RG is a 
modification of the primordial PS appearing in Eq. (3), 
as P^{k) P^{k) = P^{k)0{X - k). Inserting the 
truncated PS in Eq. (3) yields a generating functional 
Zx[Ja, Kb] P^] = Z[Ja, Ki,] P^], that describes a ficti- 
tious Universe, whose statistics of initial data is modified 
by suppressing all fluctuations with wavenumber larger 
than A. On the other hand, the dynamical content, en- 
coded in the linear propagator and in the structure of the 
interaction, is left unchanged. 

In the A ^ 00 limit, all the fluctuations are included, 
and we recover the physical situation. Increasing the 



3 



d/dA = 1/2 . 

a b i 



...9 



/ 

a b 



d/dA 



FIG. 2: RG equation for the propagator Gah,\ 
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FIG. 3: RG equation for $a6,A 



cutoff from A = to A ^ oo, tfie linear and non-linear 
effect of fluctuations of higher and higher wavenumber is 
gradually taken into account. This process is described 
by a RG equation [5, 6] which can be obtained by taking 
the A derivative of 



dxZx 



^ j driadribd^c[5{X-q)P°^{q)S{ria)S{rib)- 



5Kb5Ka 
(5) 

Taking successive derivatives w.r.t. the sources of this 
master equation - or of the analogous ones for W and 
F - one can obtain the RG evolution of any physical 
quantity. The structure of the equation is such that the 
evolution of a correlation function of order n involves all 
correlations up to order n + 2. The RG equation for the 
propagator and for the function ^ab,x are represented in 
Figs. 2 and 3, respectively, where the thick lines indicate 
full propagators, the dark box is the full PS Pab,x^ dark 
blobs are full IPX functions, and the crossed box is the 
RG kernel obtained by deriving w.r.t. A the function 
multiplying the PS in P^, that is 

(6) 

A recipe can be given to obtain the RG equation for 
any given quantity [6]: 

i) write down the 1-loop expression for the quantity of 
interest, obtained using any needed vertex^ (for instance, 
in Figs. 2, 3 we have not only the vertex but also 

XX^^ XW^ XXWi which vanish at tree- level); 

a) promote the linear propagator, the PS and the ver- 
tices appearing in that expression to full, A-dependent 
ones; 

in) take the A-derivative of the full expression, by 
considering only the explicit A-dependence of the step- 
function contained in P/. 

It should be emphasized that the RG equations ob- 
tained following these rules are exacts in the sense that 
they encode all the dynamical and statistical content of 
the path-integral (3) or, equivalently, of the continuity 



and Euler equations supplemented by the initial PS. Dif- 
ferent approximation schemes can be attempted, includ- 
ing of course PT, which can be recovered by using quanti- 
ties up to the /-th loop order in the RHS to get the / + l-th 
one by performing the A-integration. However, our RG 
equations are most indicated for non-perturbative resum- 
mations. 

As a first step, we note that the full 3 and 4-point 
functions appearing in the RG equations for Gah,\ and 
^ah,x are also A-dependent quantities, which can be com- 
puted by RG equations, also derived from Eq. (5). These 
equations depends, in turn, on full, connected, and A- 
dependent functions up to 5 (for the 3-point function) or 
6 (for 4-point ones) external legs, which also evolve ac- 
cording to RG equations. Approximations to the full RG 
flow then amount to truncating the full hierarchy of cou- 
pled differential equations, and using some ansatz for the 
full n-point functions appearing in the surviving equa- 
tions. We will approximate the full RG flow by keeping 
the running of the 2-point functions (propagator and PS) 
and keeping the tree-level expression for the trilinear ver- 
tex XW' ^his approximation, only the first diagrams 
on the RHS of Figs. 2 and 3 contribute to the running. 

The 1-loop result for the propagator [2] can be recov- 
ered by using the tree-level expressions for the kernel 
Kah,\ and for the propagators on the RHS of Fig. 2. 

Using running propagators on the RHS of Fig. 2, while 
keeping the kernel at tree- level, we get a RG equation 
which can be analytically integrated in the k ^ q = X 
limit, 

dxGab,x{k'i^a,r]b) = -Gab,x{k'^Va,r]b) 



X J d^qS{X-q) 



P\q) 



(7) 



having used the property ix/7e/^(— k, q, k — q) — 
5eg{k/2q) cosk-q, valid in this limit. Imposing the initial 
condition Gab,x=o{k, r]a, r]b) = gabiVa, Vb), and integrating 
up to A = oo, one gets 



Gab{k, r]a, Vb) =gab{Va^ Vb)e~ 



(e^a_e^b)2 



(8) 



where is the velocity dispersion, defined as = 
(1/3) / d^ c\P^ {q) / q^ . Two comments are in order. The 
RG improvement discussed here has a clear interpreta- 
tion in terms of PT. Indeed, as shown in [3], the result 
in Eq. (8) can be obtained also by resumming an infinite 
class of diagrams. It is amazing how the same result, that 
in PT requires a careful control of the combinatorics, is 
here obtained by a simple, 1-loop, integration. The sec- 
ond comment has to do with the dramatic modification 
of the UV behavior of the resummed propagator w.r.t. 
the linear one, and on its impact on the RG flow. In- 
deed, when the propagator of Eq. (8) is employed in the 
kernel Kab,x^ an intrinsic UV cut-off is provided to the 
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RG flow: fluctuations with large momenta are exponen- 
tially damped, so that the RG evolution freezes out for 
A /cFy. This is a genuinely non-perturbative ef- 

fect, which is masked if one considers PT at any flnite 
order. In other words, the full equations of motion are 
much better behaved in the UV than their perturbative 
approximations . 

Improving further on the approximations leading to 
Eq. (8), we relax the k ^ \ condition, keeping the full 
momentum dependence of the tree-level vertex, and us- 
ing A-dependent propagators also in the kernel (6). For 
the full propagator we use the ansatz Gae,x{k^ si, 52) = 
Hae,x{k) exp[— /c^cr^ {e ^-e which we have fac- 

tored out the leading time behavior of Eq. (8). Inserting 
the expression above in the RG equations we obtain a 
closed system of four coupled differential equations, one 
for each component of Hae, x- In order to simplify it fur- 
ther, we consider the two combinations Gai,x + Ga2,x-, 
{a = 1, 2), therefore, for each external momentum /c, we 
solve a system of two coupled differential equations. Full 
details will be given in [6]. The exponential damping of 
Eq. (8) is exhibited also by our more reflned scheme. 

The RG evolution of the propagator governs that of the 
-^ab X contribution to the full PS. On the other hand, the 
evolution of the other contribution, P^^ com- 
puted by solving the equation in Fig. 3. We approximate 
the exact RG equations along the same lines we follow in 
the computation of the propagator. In particular, the 77 
integrations in the expression for P^^ require an ansatz 
for its 'time' dependence. In analogy with Eq. (8), we use, 
again on the RHS of the RG equation, P^l ^{^^ Va^ Vb) = 

compute the PS at equal 'times', rja — = 7^, with the 
initial condition P^^ x=oi^) — ^ • We consider a spatially 
flat ACDM model with n% = 0.7, = 0.046, h = 0.72, 
Us = 1. The primordial PS P^ is taken from the out- 
put of linear theory at Zin = 35, as given by the CAMB 
Boltzmann code [7]. 

In Fig. 4 we plot our results for the PS (solid lines), in 
the momentum range relevant for the BAO, at 2: = 0, 1, 
and 2. The short-dashed lines correspond to the linear 
theory and the long-dashed ones to 1-loop PT (which, at 
z = has been truncated for k ^ 0.17 h/Mpc, where P^ 
takes negative values, signaling the breakdown of the per- 
turbative expansion). The black squares are taken from 
the numerical simulations of Ref. [8]. To enhance the 
BAO feature, each PS has been divided by the linear one, 
in a model without baryons [9]. In the peak region, our 
RG results agree with those of N-body simulations to a 
few percent accuracy down to redshift z = 0, where linear 
and 1-loop perturbation theory badly fail. Thus, the dy- 
namical behavior in this momentum range appears to be 
captured fairly well by the approximations implemented 
by our approach, namely the 'single stream approxima- 



tion', leading to Eq. (1), and the non-linear corrections of 
z=2 z=l 




FIG. 4: The power-spectrum at 2; = 2, 1, 0, as given by the 
RG (solid line), linear theory (short-dashed), 1-loop PT (long- 
dashed), and the N-body simulations of [8] (squares). 



the two-point functions only, i.e. the propagator and the 
PS. The RG performance can be systematically improved 
by increasing the level of truncation of the full tower of 
differential equations, the next step being the inclusion 
of the running of the trilinear vertex. 
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